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We present some results of a numerical investigation of the 't Hooft model with 2, 3 and 4 colors, regularized 
on the lattice with overlap fermions. 



Several years ago, in a pioneering investigation, 
't Hooft studied U{Nc) theories in the limit Nc — > 
oo with Ct = g'^Nc kept constant [|]. In two di- 
mensions this led him to a model of QCD that ex- 
hibited important features of the full theory. Re- 
cently, the overlap formulation has made it possi- 
ble to introduce a lattice Dirac operator D which 
in the massless limit preserves a lattice form of 
chiral symmetry at finite lattice spacing |^ . This 
development, together with a renewed interest in 
't Hooft 's results, prompted us to perform a nu- 
merical investigation of a class of two-dimensional 
non-Abelian models with overlap lattice femions. 
Precisely, we simulated models of QCD2 with 
SU{Nc) and U{Nc) color groups for Nc = 2,3,4, 
in the 't Hooft limit Ct = g^Nc = const. Since 
in two dimensions only the U{Nc) models exhibit 
a U{1) axial anomaly, we studied both classes of 
models in order to analyze the behavior of the 
overlap operator in presence and absence of the 
axial anomaly. Our systems are small enough 
that it is possible to compute exactly in the over- 
lap formulation the propagator and det{D), 
following the scheme used in We computed 
the masses of the pions as a function of the quark 
mass TO and, for both classes of models, we found 
the expected functional dependence on to and 
universality in Nc- In the U{Nc) case the me- 
son spectrum exhibits also a flavor singlet parti- 
cle (77'), whose mass is known analytically, and 
which we computed by following the approxi- 
mate method of Ref. [|j and by using the Witten- 
Veneziano relation. The details of our simulation. 



*Presented by F. Berruto. Work supported in part by 
grant DE-FG02-91ER40676 and by INFN. 



as well as additional results that we could not in- 
clude here because of lack of space, will be pre- 
sented in a separate article Q . 

1. The 't Hooft Model with Overlap 
Fermions 

The 't Hooft model action in the overlap regu- 
larization reads 
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where [/^^ is the usual Wilson plaquette, a is the 
lattice spacing and are Nf fermion fields with 
mass TO and (3 = 2Nc/{agf. D ^ {l + V)/a is 
the massless Neuberger-Dirac operator, V being 
the unitary component of Dw — where Dw 
is the usual Wilson-Dirac operator. 

For TO = the action (H) is invariant under a 
continuous symmetry Sij) = 75(1 — aD)il}, 5tJj = 
■075, i.e. a form of chiral symmetry that holds at 
finite lattice spacing. The Nf = 2 models exhibit 
an isotriplet of particles (tt), whose masses vary 
with the quark mass as ||l[] 1^ 
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U{Nc) models (2) 
+ ... , SU{Nc) models (3) 



Eq. (j|) has been derived for Nc 00 while Eq. (j|) 
also involves a semiclassical WKB approximation. 

In the U{Nc) models in the chiral limit M^y — 
Nf g'^ /tt. The Witten-Veneziano formula for the 
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U{Nc) models gives 
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2. The meson spectrum 



We performed extensive simulations with the 
exact Neuberger-Dirac operator on a 18 x 18 lat- 
tice at Ct = 4/3, i.e. /3 = 6, 13.5, 24 for 
Nc = 2,3,4 respectively. We computed for 
TO = 0.04,0.05,0.06,0.07,0.08,0.1 so that the re- 
lation NxMtt > 4 holds. The effects of dynam- 
ical fermions have been included by weighting 
the observables with the appropriate power of the 
fermion determinant. The smallness of the lattice 
warrants this procedure (cfr. j^), which would 
lead to unacceptable variance on larger systems. 

We computed the meson propagators projected 
over zero momentum from the fermion propaga- 
tors and extracted the meson masses in a stan- 
dard manner. The errors have been estimated 
with the jacknife analysis. 

In table (|^) we report the pion masses for the 
U{Nc) and SU (Nc) models with two flavors of dy- 
namical fermions and in the quenched approxima- 
tion. The quenched results get closer to the un- 
quenched two flavors results when Nc gets larger, 
as one would expect since the large Nc limit cor- 
responds to discarding the internal quark loops. 
For the U{Nc) case Ec^.f^) motivated us to do the 
fit 
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We obtained A ^ -0.012(17), B = 1.57(19) for 
Nc = 2, A^ -0.006(4), B = 1.53(4) for iV^ = 3 
and A ^ -0.006(8), B = 1.66(8) for Nc = 4. 

For the SU{Nc) models, following Eq.(||), we 
performed the fit 
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We obtained A = 0.022(26), B = 1.524(120) for 
Nc^2, A^ 0.045(43), B = 1.548(166) for Nc = 
3 and A = -0.009(40), B = 1.788(11) for Nc = 
4. In the SU (Nc) case we performed also the fit 
I g\fWc = C{m/ g\/WcY and we obtained 7 
compatible with 2/3 for Nc = 2, 3, 4. Figures(0,|) 



U{Nc) models 




Nc^2 


Nc^i 


TVc = 4 


rriq/g^/Nc 


APjg^Nc 


Ml/g^Nc 


APjg^Nc 




Nf = 


2 




0.0346 


0.045(10) 


0.047(2) 


0.055(6) 


0.0433 


0.056(9) 


0.059(2) 


0.067(5) 


0.0519 


0.068(8) 


0.072(3) 


0.079(4) 


0.0606 


0.080(7) 


0.085(3) 


0.093(3) 


0.0692 


0.094(5) 


0.098(3) 


0.107(2) 


0.0866 


0.124(4) 


0.128(3) 


0.138(1) 


quenched {Nf — 0) 


0.0346 


0.069(1) 


0.063(1) 


0.057(1) 


0.0433 


0.083(1) 


0.076(1) 


0.071(1) 


0.0519 


0.098(1) 


0.091(1) 


0.085(1) 


0.0606 


0.113(1) 


0.106(1) 


0.100(1) 


0.0692 


0.128(1) 


0.121(1) 


0.116(1) 


0.0866 


0.161(1) 


0.153(1) 


0.148(1) 
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iV, = 2 ^ 


iV, = 3 ^ 


Nc^^ ^ 


(TO,/g7v|)2/3 


M^/gNl 


M^/gNc^ 


M^/gNl 




Nf = 
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0.1062 


0.183(14) 


0.206(23) 


0.181(28) 


0.1233 


0.211(12) 


0.237(23) 


0.211(26) 


0.1392 


0.236(10) 


0.263(21) 


0.239(24) 


0.1543 


0.258(8) 


0.286(19) 


0.266(22) 


0.1686 


0.279(7) 


0.307(16) 


0.292(21) 


0.1957 


0.321(5) 


0.348(10) 


0.341(17) 


quenched {Nf — 0) 


0.1062 


0.200(2) 


0.212(1) 


0.217(1) 


0.1233 


0.227(2) 


0.241(1) 


0.246(1) 


0.1392 


0.252(2) 


0.267(1) 


0.273(1) 


0.1543 


0.276(1) 


0.292(1) 


0.298(1) 


0.1686 


0.299(1) 


0.316(1) 


0.322(1) 


0.1957 


0.343(1) 


0.361(1) 


0.367(1) 



Table 1 

give evidence of universality in Nc for the pion 
masses. 

In the U[Nc) models we computed also M^/. 
The rj' propagator is given by the difference be- 
tween the connected and disconnected terms in 
the correlator and, because of cancellations, the 
errors are larger than in the tt case. There- 
fore, we resorted the method proposed in Q 
to compute the 77' mass. The method exploits 
the quenched two-loop disconnected Fg ^'^'-'P and 
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Figure 1. M^/g'^N^ vs. nig/g^/l^ for Nf = 2, 
U{Nc) models 
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Figure 2. M^/gy/Nl vs. {mq/ gy/W^f/^ for Nf 
2, SU{Nc) models 
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tions to the 77' propagator: 
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As one can see from table (|^) M^, computed us- 
ing Eq. and N f ~ 1 compares remarkably well 
with the analytic result. 

An alternative calculation of M^/ can be done 
by using the Witten-Veneziano formula. In 
the quenched approximation we computed = 
0.774(3), 0.956(3), 1.115(5) for = 2,3,4 re- 
spectively. Since in our calculations we diag- 
onalize the full Neuberger-Dirac operator, we 
could determine the number of zero modes, and 
thus the topological charge, for each individual 
gauge configuration. We thus obtained (Q^) /V = 
0.0258(16), 0.0298(23), 0.0319(37) for = 
2,3,4 respectively. (Notice that with increasing 
Nc these values get closer and closer to the ana- 
lytic value {Q'^)/V = Ct/^i^'^ = 0.0337. Building 
on the trade-off between spatial and internal de- 
grees of freedom that prompted the introduction 
of the Nc ^ 00 single plaquette model Q, one 
could heuristically argue that the theory should 
be less affected by finite size effects increasing 
Nc. A rigorous justification of this argument, 
however, would require a detailed analysis of dis- 
cretization and finite size effects that goes beyond 
the scope of our investigation.) From /^r and 
{Q^)/V, we use Eq.(E|) to calculate Af^/. The re- 
sults, also in table (g), are consistent for iVc = 3 
and 4 with those obtained with the method of 
Ref. M and with the analytical formula. 



MlJg^N, 


iVc = 3 


A^c = 4 


Analytic 


0.159 


0.106 


0.079 


Numerical from Eq.( 




0.159(8) 


0.111(7) 


0.084(8) 


Numerical from Eq.(y) 


0.129(6) 


0.097(7) 


0.076(6) 



Table 2 



quenched one-loop connected t\ ^'^'^P contribu- 



REFERENCES 

1. G. 't Hooft, Nucl. Phys. B 75, 461 (1974). 

2. H. Neuberger, Phys. Lett. B 417, 141 (1998); 
Phys. Lett. B 427, 353 (1998). 

3. L. Giusti, C. Hoelbling and C. Rebbi, Phys. 
Rev. D 64, 054501 (2001). 

4. H. W. Hamber, E. Marinari, G. Parisi and 
C. Rebbi, Nucl. Phys. B 225, 475 (1983). 

5. F. Berruto, L. Giusti, C. Hoelbling and 
C. Rebbi, in preparation. 

6. C. J. Hamer, Nucl. Phys. B 195, 503 (1982). 

7. T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 
1063 (1982). 



